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The well-known results concerning a dilute Bose gas with the short-range repulsive interaction 
should be reconsidered due to a thermodynamic inconsistency of the method being basic to much of 
the present understanding of this subject and nonrelevant behaviour of the pair distribution function 
at small boson separations. The aim of our paper is to propose a new way of treating the dilute Bose 
gas with an arbitrary strong interaction. Using the reduced density matrix of the second order and 
a variational procedure, this way allows us to escape the inconsistency mentioned and operate with 
singular potentials like the Lennard- Jones one. All the consideration concerns the zero temperature. 



I. INTRODUCTION AND BASIC EQUATIONS 

It is well-known that to investigate a dilute Bose 
gas of particles with an arbitrary strong repulsion (the 
strong-coupling regime), one should go beyond the Bo- 
goliubov approach Q (the weak-coupling case) and treat 
the short-range boson correlations in a more accurate 
way. An ordinary manner of doing so is the use of the 
Bogoliubov model with the "dressed" , or effective, inter- 
action potential containing "information" on the short- 
range boson correlations (see Ref. |^). Below it is demon- 
strated that this manner leads to a loss of the thermody- 
namic consistency. To overcome this trouble, we propose 
a new way of investigating the strong-coupling regime 
which concerns the reduced density matrix of the sec- 
ond order (the 2-matrix) and is based on the variational 
method. 

The 2-matrix for the many-body system of spinless 
bosons can be represented as H: 



P2(ri,r^;ri,r2) 



F2(ri,r2;r;,r'2) 



(1) 



N{N-l) ' 

where the pair correlation function is given by 

F2(ri,r2;r;,4) = {i;\v^)i,\v2mv'^)i^{v[)) . (2) 

Here ?/'(r) and V'^(r) denote the boson field operators. 
Recently it has been found that for the uniform sys- 
tem with a small depletion of the zero-momentum state 
the correlation function (^) can be written in the ther- 
modynamic limit as follows: 

i^2(ri,r2;rl,r'2)=n2y,*(r)^(r') 

+2*^0 j -fl^nq (p;/2W<^q/2(r') 

X exp{iq(R' - R)}, (3) 

where r = ri — r2, R = (ri -I- r2)/2 and similar rela- 
tions take place for r' and R', respectively. In Eq. ^ 
uq = Nq/V is the density of the particles in the zero- 
momentum state, Uq = (aJjCq) stands for the distribution 
of the noncondensed bosons over momenta. Besides, (p{r) 



is the wave function of a pair of particles being both con- 
densed. In turn, 'y3q/2(r) denotes the wave function of the 
relative motion in a pair of bosons with the total momen- 
tum /iq, this pair including one condensed and one non- 
condensed particles. So, Eq. (^) takes into account the 
condensate-condensate and supracondensate-condensate 
pair states and is related to the situation of a small de- 
pletion of the zero-momentum one-boson state. For the 
wave functions ip{r) and <Pp(r) we have 

(p{r) = l+'ijj{r), (^p(r) = V2cos(pr) + V'p(r) (p 7^ 0) (4) 

with the boundary conditions ip{r) — s- and i/jpi^) — > 
for r oo. The functions ip{r) and ^/'p(r) can explicitly 
be expressed in terms of the Bose operators and ap 



ip{k) = / VC?-) exp(-ikr) d^r = (aka_k)/n-o, (5) 



V'p(k) = / ^p{r) exp{~ikr) d-^r 



2nn 



«2p 



(6) 



Having in our disposal the distribution function rik and 
the set of the pair wave functions f{r) and </?p(r), we are 
able to calculate the main thermodynamic quantities of 
the system of interest. In particular, the mean energy 
per particle is expressed in terms of Uk and g{r) via the 
well-known formula 



(27r)3 n 2 

where Tk = m is the one-particle kinetic energy, 

n = N/V stands for the boson density and the relation 

ff(r) =i^2(ri,r2;ri,r2)/n2. (8) 
is valid for the pair distribution function g{r). 
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II. THE BOGOLIUBOV MODEL 

The starting point of our investigation is the weak- 
couphng regime which imphes weak spatial correlations 
of particles and, thus, is characterized by the set of the 
inequalities 



\^{r)\ « 1, |^p(r)| « 1 



(9) 



Specifically, the Bogoliubov model corresponds to the 
choice |,| 



|^(r)|«l, Vp(r)-0 



(10) 



Besides, owing to a small depletion of the Bose conden- 
sate (n — no)/n we have for the one-particle density ma- 
trix^^i(r) = (^t(i.i)V,(r2)): 



Fi(r) 



nk 



(27r)3 n 



exp(ikr) 



< 



no 



< 1. 



So, investigating the Bose gas within the Bogoliubov 
scheme, we have two small quantities: and Fi{r)/n. 
This enables us to write Eq. (||) with the help of (|]) as 
follows: 



n 



(2^ 



nfcexp(ikr), (11) 



where we restricted ourselves to the terms linear in V'l'') 
and Fi{r)/n and put '>p*{r) — ^{r) because the pair wave 
functions can be chosen as real quantities. Equations for 
^(k) and rik can be found varying the mean energy (|^) 
with Eq. ( pi] ) taken into account. However, previously 
one should realize an important point, namely: Uk and 
^{k) can not be independent variables. Indeed, when 
there is no interaction between particles, there are no spa- 
tial particle correlations either. So, V'(^) = and, since 
the zero-temperature case is considered, all the bosons 
are condensed, Uk = 0. While "switching on" the inter- 
action results in appearing the spatial correlations and 
condensate depletion: '4'{k) ^ together with Uk ^ 0. In 
the framework of the Bogoliubov scheme is related 
to Uk by the expression 



nk{nk + 1) = nltj?{k) 



(12) 



Indeed, the canonical Bogoliubov transformation jl] im- 
plies that 

Ok = Wfettk + Wfcalk' aj^ = Ufcaj^ -I- ?;fca_k, (13) 
where 



(14) 



At zero temperature (aj^ak) — and, using Eqs. (|^) and 
(O) we arrive at 



nk=vl, ij^ik) = UkVk/no. 



(15) 

With Eqs. (|lj) and (|l|) one can readily obtain Eq. (|l|). 

Now, let us show that all the results on the thermo- 
dynamics of a weak-coupling Bose gas can be derived for 
the Bogoliubov scheme with variation of the mean en- 
ergy under the conditions ( [TTI ) and (12). Inserting 
Eq. ( pll ) into Eq. (|^) and, then, varying the obtained 
expression, we arrive at 



5e = 



(fk 



Tk + n$(fc) ) — + nMk)5i>{k) 



Relation (|12|) connecting ipik) with Uk results in 
{2nk + l)dnk , V^(fc) f d^q 



Sipik) 

2nl%l)(k) 

where the equality 

n — no + 



no 



(27r) 



:^Snq, 



d^k 
(2^ 



nk 



(16) 



(17) 



(18) 



is taken into consideration. Setting fe = and using 
Eqs. ( p6| ) and (|l^), we derive the following expression: 

^2 



2Tki;{k) 



$(fc) 



2nk) + 2nijj{k) 
d^q 



n 
no 



{2^y 



<i>((j)^(g) . (19) 



Here one should realize that Eq. ( |19| ) is able to yield re- 
sults being accurate only to the leading order in (n — 
no)/n because the used expression for g{r) given by Eq. 
( pT| ) is valid to the next-to-leading order So, Eq. ( p^ ) 
should be rewritten as 

-2Tki^{k) = $(fc)(l + 2nk) + 2n^{k)<^{k). (20) 

Equation ( pO| ) is an equation of the Bethe-Goldstone type 
or, in other words, the in-medium Schrodinger equation 
for the pair wave function. As 2$(fe)(nfe -I- n%l}{k)) is the 
product of the Fourier transforms of $(r) and n{g(r) — l), 
we can rewrite Eq. ( pO| ) in the more customary form 

^VV(^) = Hr) +nj mv - y|)(<?(y) - l)d^y. (21) 

The structure of Eq. (^) is discussed in the papers [S,^ . 
Here we only remark that the right-hand side (r.h.sj of 
Eq. ( |2l] ) is the in-medium potential of the boson-boson 
interaction in the weak-coupling approximation. The sys- 
tem of equations ( p^ ) and ( [20|)_c an easily be solved, which 
leads to the familiar rcsultsfQ: 



nk = -z 




Tk + n^jk) 

n + 2nTk^{k) 
$(fc) 



- 1 



+ 2nTk<^{k) 



(22) 
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III. A DILUTE BOSE GAS WITHIN THE 
BOGOLIUBOV MODEL 



The derived integral is readily calculated. The result is 
given by 



As it was mentioned, the aim of our paper is investi- 
gation of the case of a dilute Bose gas with an arbitrary 
strong repulsion between bosons. So, considering a di- 
lute Bose gas in the weak-coupling approximation can be 
a good exercise providing us with useful information. Let 
us investigate the thermodynamics of a dilute Bose gas 
within the Bogohubov model. With Eqs. (|^), and 
(E2) we derive 



n~,^. 1 f (Pk 



T2 + 2nTfe$(fc)-rfc-n$(fc) 



(23) 



The well-known argument of Landau (see the footnote 
in Ref. j^] and discussion in Ref. ||^) testifies that the 
properties of dilute quantum gases are ruled by the scat- 
tering length. Within the Bogoliubov model this length 
is usually assumed to be equal to m$(0)/47r?i^. If so, 
when expanding e in powers of the boson density n, one 
could replace $(fc) by $(0) in Eq. (p3[), introducing the 
low-momentum approximation. However, this leads to 
a divergency because at large k the integrand behaves 
as — n^$^(A;)/2Tfc. To properly calculate the integral in 



Eq. (|23|), we should rewrite Eq. 
form: 



|23|) in the following 



e = 



(fk 

(2^)3 2Tk 



(24) 



where 



1 f „ 47rfc2 
2n J (27r)3 





(v^ 



x[^jT^ + 2nTk^k) -Tk- n$(fc) 



Now, substituting k — (2mnyy/'^/h in the integral, we 
obtain the expression 



V2 
47r2 



3/2 



oo 

J y + 2$(^/2mny/n) 



-,3/^-$(V2;;^A),v. + £!(v^^| 

which at sufficiently small n may be rewritten as 

3/2 



y/2 / mn \ ' 



X /■ dy{y^y + 2$(0) - y^^ - $(0)yi/2 + 



2Vy 



I = 



157r2 



In turn, the first term in the r.h.s. of Eq. (| 
represented as 



(25) 



can be 




- / ^(°)(r)$(r)dV, 



(26) 



where (p'-"-' is the solution of Eq. ( pi] ) in the limit n ^ 0. 
This is nothing else but the Schrodinger equation in the 
Born approximation. According to the relations (^5|) and 
( ^ ) we has to conclude that the scattering length in the 
case of interest is expressed in the form 



as 



47r?i^ 



ip("\r)<i>{r)d^r. 



(27) 



One can easily be convinced that <&(0) can not be repre- 
sented only in terms of as and, hence, the dependence 
on the shape of the interaction potential appears in the 
series expansion for the mean energy in the first correc- 
tion to the term iTrfi^aBn/rn. To rewrite our result for e 
in a graphic form, we introduce one more characteristic 
length 6 > which obeys the relation 



h = -- 



d?k <P^{k) 



A7:h^ J (2^)3 2Tk 



(28) 



where (r) = ip^*^'>{r) — 1. Further, with the help of 
Eqs. (p3)-(p3), we arrive at 



2Trh^aBn f 128 



1 + ^) 

2aB 



(29) 



here the condition 6 <C as is of use. It is not difficult 
to see that the expression (|2^) taken with negative sign 
is the next correction to the scattering length calculated 
within the Born approximation. As to the Eq. (^^, it is 
related to the next-to-Born approximation. Stress that in 
the Bogoliubov model the energy term n<I>(0)/2 is treated 
as the major one ||l|, which implies that the condition 
6 <C as is fulfilled. This qualitative criterion can be 
written as 



$(0) > 



d?k 

(27r)3 2Tu ■ 



(30) 
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Beyond this inequahty the model may be thermodynam- 
ically unstable. In particular, the opposite case 



$(0) < 



(27r)3 2Tk 



leads to the negative scattering length which at suffi- 
ciently low densities results in —d^E/dV^ — dp/dV > 0. 

Thus, investigated within the Bogoliubov model, the 
thermodynamics of a dilute Bose gas is ruled by the 
scattering length only in the zero-density limit. While 
the next-to-leading term in the series expansion given by 
Eq. (I29I ) depends on the shape of the interaction, which 
is expressed in appearance of the additional character- 
istic length b. This conclusion differs from the results 
of papers |^ according to which the series expansion for 
e taken to the same order as that of Eq. (p9|), is fully 
determined by the scattering length. To clarify the sit- 
uation concerning this difference, we should go to the 
strong-coupling regime. 



IV. THE STRONG-COUPLING REGIME 

Now, after the detailed investigations of the Bogoli- 
ubov model within the scheme proposed, we are able to 
demonstrate that the investigation of the strong-coupling 
case based on the Bogoliubov model with the effective 
boson-boson interaction, results in a loss of the thermo- 
dynamic consistency. Indeed, as it was shown in the pre- 
vious section, any calculating scheme using the basic re- 
lations of the Bogoliubov model (pT|), ( p^ ) conclusively 
leads to Eqs. (|2C|)-(|2^) provided this scheme does yield 
the minimum of the mean energy. In this case Eqs. (pO|)- 
( p^ ) certainly includes the quantity <I>(r) which is the 
"bare" interaction potential appearing in Eq. (0). The 
use of the Bogoliubov model with the effective interac- 
tion potential substituted for <I>(r) can in no way disturb 
the relations given by Eqs. (|ll|) and (12). And Eq. (^ is 
the same in both the weak- and strong-coupling regimes. 
Thus, any attempts of replacing <I>(r) by the effective 
"dressed" potential without modifications of Eqs. ( [ll| ) 
and (|l^) results in a calculating procedure which docs 
not really provide the minimum of the mean energy. It is 
nothing else but a loss of the thermodynamic consistency. 
We remark that we do not mean, of course, that the t- 
matrix approach or the pseudopotential method can not 
be applied in the quantum scattering problem. It is only 
stated that the usual way of combining the ladder dia- 
grams with the random phase approximation faces the 
trouble mentioned above. Though our present investi- 
gation is limited to the consideration of the many-boson 
systems, the derived result gives a hint that the similar 
situation is likely to take place in the Fermi case, too. In 
this connection it is worth noting the problem associated 
with the lack of self-consistency of the standard method 
of treating the dilute Fermi gas M . 



The strong-coupling regime is characterized by sig- 
nificant spatial correlations. So, Eq. (10) resulting in 
Eq. (|l^) is not relevant for an arbitrary strong repul- 
sion between bosons at small separations when we have 
V'(O) = -1, i/'p(0) = -%/2 (see Refs. |,|). Therefore, 
to investigate the strong-coupling regime, Eq. (llj) should 
be abandoned in favor of Eq. (||). Expression (|) is accu- 
rate to the next-to- leading order in (n — uq) /n. So, using 
Eqs. (0) and (||), we can write 



g(r) =ip (r) + - 
n 



(27r)3"' 



{^l/,{v)-^\r)). (31) 



Let us now perturb ip{k) and n{k). Working to the first 
order in the perturbation and keeping in mind conditions 
(|l|) and (jsif), from Eq. (|) we derive: 



2Tfei/'(fc) = U{k){l + 2nk) + 2n^:{k)U'{k) (32) 



with 



and 



U{k) ^ / (p(r)cf>(r) exp(-ikr)d^r 



(33) 



(34) 



Using Eqs. (|33D, ( |34| ) as well as the relation V'k(r) 
V2^{r) {k 0) (see the boundary conditions (^) [||, 
we obtain U{0) U'{0). This implies that the system 
of Eqs. (|l^ ) and ( |3^ ) is not able to yield the relation 
Ilk oc 1/fc (fc ^ 0) following from the "1/fc^" theorem 
of Bogoliubov for the zero temperature Indeed, let 
us assume Uk 00 for fc ^ 0. Then, from Eq. (12) at 
n — hq we find n\^p{k)\/nf^ — > 1 when fc ^ 0. On the 
contrary, Eq. @) gives n\if{k)\/nk Lf{0)/U'{0) ^ 1 
for fc ^ 0. So, consideration of the Bose gas based on 



Eqs. (H) and (|T2|) does not produce satisfactory results. 
Nevertheless, it is worth noting that Eq. (|3^) has an im- 
portant peculiarity which differentiate it from Eq. ( [20| ) 
in an advantageous way. The point is that in both the 
limits n — > and fc — > cx3 Eq. (p3) is reduced to 



(35) 



As it is seen, this is the exact "bare" (not in-medium) 
Schrodinger equation, other than its Born approximation 
following from Eq. (|2^). Thus, we can expect the line of 
our investigation to be right. 

As it was shown in the previous paragraph, an ap- 
proach adequate for a dilute Bose gas with an arbitrary 
strong interaction can not be constructed without mod- 
ifications of Eq. (|2|). This is also in agreement with a 
consequence of the relation 



|(aka_k)P < (aka|^)(a^^^a_k) 



(36) 
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resulting from the inequality of Cauchy-Schwarz-Bogo- 
liubov Q 

\{AB)\^ < {AA^){B^B). 

With Eqs. (^) and (|6|) one can easily derive nQtp'^{k)< 
nk{nk + l)- Thus, it is reasonable to assume that Eq. (|l3) 
takes into account only the condensate-condensate chan- 
nel and ignores the supracondensate-condensate ones. 
Now the question arises how to find corrections to the 
r.h.s. of Eq. (p^. At present we have no regular proce- 
dure allowing us to do this in any order of (n — no)/n. 
However, there exists an argument which makes it possi- 
ble to realize the first step in this direction. The matter 
is that the alterations needed have to produce the equa- 
tion for V-'p(k) which is reduced to the equation for tpik) 
in the limit p 0. Though this requirement does not 
uniquely determine the corrections to Eq. jl^), it turns 
out to be significantly restrictive. In particular, even the 
simplest variant of correcting Eq. ( p^ in this way, leads 
to promising results. Indeed, this variant is specified by 
the expression 



nkirik + 1) = nl ^^{k) + 2no 



703 (37) 



Eq. ( p7| ) is valid to the next-to-leading order in (n— no)/n. 
So, we may rewrite it as 

rikiuk + 1) = n^i^\k) + 2n[ .^ng(fl/^{k) - ^^^(fc)). 



(38) 



Perturbing tpW s-nd ^fc bearing in mind conditions 
(|l|) and (11), Eq. (0) gives Eq. (^) again. However, 
now U'{k) obeys the new relation 



(39) 



U'ik)^ j (ipl/^{v)~^\T))^r)cpT 

f ^(g)(V^g/2(q)-^'(g)) 



(27r)3 



V'(g) 



which significantly differs from Eq. (34). Indeed, the 
choice of the pair wave functions as real quantities im- 
plies that operating with integrands in Eqs. ( [33[ ) and 
(^) , one can exploit ■(/'p (r) — V^V' ('') small p . 

For A; ^ this provides U'{k) - U{k) =Jk = ck'^ -\ . 

Note that for fc — > oo we have tk — > —U{Q). Similar to 
Eq. (20), Eq. ( |3^ ) can yields results correct only to the 
leading order in (n — no)/n. So, it has to be solved to- 



gether with Eq. (yD where should be replaced by n 



rather than with Eq. 
relations: 



This leads to the following 



1 



nk 




Tk + nU{k) 



Tl -f 2nTkU{k) 



(40) 



U{k) 



Tl + 2nTkU{k) 



(41) 



1 at fc — > CX3. 



where Tk = Tk+ntk, with the limit Tk/Tk 

For fc ^ Eq. (|^ gives Uk ~ {\JnmU{0)/nk - l)/2, 
which is fully consistent with the "1/fc^" theorem of Bo- 
goliubov for the zero temperature |TG]. 

As it is seen, the strong-coupling regime is more com- 
plicated than the Bogoliubov one because we do not know 
the quantity U (fc) ab initio. To find it, one should solve 
Eqs. ( ^3| ) and ( pl| ) in a self-consistent manner. Equations 
( ^3|) and ( |4l| ) lead to one more interesting relation 



U{k) = $(fc) - 



1 



d^q $(|k-q|)L/(g) 
(27r)\/y2 + 2nf,C/(g) 



(42) 



which can be called the in-medium Lippmann-Schwinger 
equation for the scattering amplitude. To obtain the ex- 
pansion for the energy at low densities, we must solve 
Eq. (p^ ) at n ^ 0. Let us rewrite it in the form 

rui.\ 1 [ "^'l S(|k-q|)t/(g) 

u{k) = m-^jj^ — 1„ 

where for Ii we have 

^\ r d\ \ $(|k-q|)L^(g) _ $(|k-q|)L/(g) | 
' 2^ (2-)^\yT,^ + 2nW,) J- 

Operating with Ii in the same manner as we dealt with 
/ in the section HI and taking into account that tk ~ 
at fc = 0, for n — !■ we derive 



(43) 



From Eqs. ( [42| ) and (43) it now follows that 

d'q $(|k-q|: 



U{k)-U^°\k)=a^{k) 



{2^f 



2Ta 



x(u{q)-U^°Hq) 



(44) 



Here U'-"\k) = J Lp'-"\r)^{r) exp{-i'kr)d^r but now 
(^(0) obeys Eq. ( ^ ) rather than Eq. ( pl| ) taken in the 
limit n — > like in the section |lll[ Let us introduce the 
new quantity ^{q) = -{U{q) - U^^\q))/2Tq. Then, for 
its Fourier transform ^(r) we obtain 



— [a + C(r)) + $(r) {a + ^(r)) = 0, 



(45) 



here ^(r) — > when r — > oo. Comparing Eq. ( pq ) with 
Eq. (^5|), we find ^(r) — aifj'^^^r). Hence, for n — > we 
have 
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C/(fc) ~ J7(°)(fc) l+7(fc, 



(46) 



Here j(k,n) 1 when n 0. Besides, the relation 
[/*^°^(0) = 'i-Kh^a/m is used in Eq. (^6|), where a is the 
scattering length. 

Having in our disposal Eq. (|4^), we are able to calcu- 
late the expansion in powers of n for the condensate de- 
pletion and energy of a dilute Bose gas with an arbitrary 
strong interparticle potential. Considering the conden- 
sate depiction {n — nn)/n — 1/{2tt)^ J^°° dkiTrk^nk/n, 
with the help of Eq. ( [40|) we obtain 



no 



3^/^ 



(47) 



Notice that according to Eq. ( |46D one can expect that 
among the omitted terms in Eq. ( p7| ) there is one pro- 
portional to na^ . 

The most simple way of deriving the expansion for the 
mean energy per particle is based on using the chemical 
potential which, in the presence of the Bose condensate, 
is given by 



dV$(|r - r'|)(V/(r')V'(r')V'(r))- (48) 



This formula follows from the well-known expression 
(551 = {5 (ji — IJ,N^), where is an infinitesimal change 
of the grand canonical potential, and relation (see Ref. 
0) 



dNo 



0. 



Using the specific expressions for the scattering parts 
of the condensate-condensate and supracondensate- 
condensate pair wave functions M given by Eqs. (|) and 
(0) one can represent Eq. (Esl) in the following form: 



= noU{0) + \/2 



3 Ug CVq/ 



t/q/2(q/2). 



(49) 



here 



C^p(k) 



(^p(r)$(r) exp(— ikr)(i^r. 



Now, for n ^ (see the procedure of calculating the 
integral / in the section || one can rewrite Eq. ( ^9| ) as 



H = nU{0) 1 + 



n — hq 



n 



(50) 



Inserting Eqs. (M) and IWA) into Eq. (BQ), we arrive at 



Anh an 
fi = I 1 



(51) 



This result for the chemical potential implies, due to the 
basic thermodynamic formula /i — d{en) / dn, the follow- 
ing expansion for the mean energy per particle: 



128 



(52) 



As it is seen, the relation ( |52| ) coincides with the well- 
known result of the approach ||] being reduced to the 
Bogoliubov model with the "dressed" interaction. It 
is not a surprise because according to the conclusions 
of the section |l[ we know that the numerical factor 
128/(15y^) appears in the series expansion for the mean 
energy per particle within the Bogoliubov model (see Eq. 
(p9|)). Replacing the bare interaction potential by the 
"dressed" one results in replacing the scattering length 
as in Eq.(p9|) by its exact value a. The only problem 
of doing so concerns the parameter b. Indeed, it follows 
from Eq. (^ that substituting the hard-sphere poten- 
tial f/(0) = iTTfL^a/m for $(fc) leads to the familiar di- 
vergency (see, e.g. Ref. p. 314). This obstacle has 
been overcome with the help of the well-known argument 
of Landau (see the footnote in the paper [||) stating that 
the thermodynamics of dilute quantum gases is only ruled 
by the vacuum scattering amplitude. According to this 
reasoning one can expect that dependence on the shape 
of the interaction potential should not appear in the first 
orders of the density series expansion of the thermody- 
namic quantities. So, various regularizing procedures, 
more or less speculative, have been worked out in order 
to exclude this divergence (together with the parame- 
ter b). On the contrary, there are no problems like this 
within the approach of the present paper. Here Eq. ( [5^ ) 
is derived on the solid theoretical basis rather than with 
the help of Landau's argument. In spite of its reasonable 
character, it needed to be corroborated, and the results of 
this paper given by Eqs. (^), (^) and (|5^) have proved 
the validity of Landau's argument beyond any inconsis- 
tencies and divergencies. In the weak-coupling case when 
|$(r)| <C 1, the energy per particle calculated within our 
scheme is expressed by Eq. ( ^2| ) with a replaced by a^. 
So, the appearance of the parameter b in the results of 
the section |l| is an artifact following from the neglect of 
scattering in the supracondensate-condensate pair wave 
channel. 

The divergence mentioned in the previous paragraph is 
not typical of the strong-coupling perturbation theory for 
the many-boson systems but results from, say, the weak- 
coupling spirit of the approach of Ref. . A simple way 
to be convinced of this is to consider the spatial boson 
correlations. Taken to the lowest-order with respect to 
the density, the structural factor (see the last paper in 
Ref. D) is of the form 



S{k) 



Tk 



T^ + 2nTkU("){k) 



(53) 



By definition we have 
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g{r) = 1 + 



-{S{k) - l)exp(ikr). 



(54) 



n J (27r)3 

Using Eqs. (m) and (p4), for n ^ one can readily find 



5(r)^l + 27/;(")(r), 



(55) 



where ^^^'>{r) obeys Eq. (Bq). This result answers the ap- 
proximation ( pi] ) while ■0'^(r) is not related to the weak- 
coupling regime and obeys the exact "bare" Schrodinger 
equation. In the situation $(r) ^ oo for r — s- one has 
■0(0) (7. — 0) —1, which implies, according to Eq. (55), 
g{r = 0) ^ —1 for n ^ 0. It is not consistent with the 
physical sense of g{r) and has nothing to do with the 
strong-coupling case corresponding to Eq. (^) when for 



g(r)-. (l + ^(°)(r))' 



Notice that the zero-density limits for the thermody- 
namic quantities of a strongly interacting dilute Bose gas 
were first found in the Bogoliubov original paper 0| : 

{n - no)/n ^ 0, <?(r) ^ ((p(°)(0)', e/n - C/(°)(0)/2. 

At last, we remark that due to the incorrect picture of 
the spatial boson correlations found in papers , one can 
expect significant alterations for the spectrum of the ele- 
mentary excitations too. However, to clarify these correc- 
tions we should conclusively solve the problem concerning 
relation between the momentum distribution and scatter- 
ing parts of the pair wave functions. Indeed, it has been 
mentioned that there exist various possibilities of gener- 
alizing Eq. (|l^ ) so as to obtain the equation for V'p(k) 
being reduced to the equation for ip{k) in the limit p — > 0. 
These possibilities result in the same series expansions for 
the thermodynamic quantities (^, (|5l|) and (|2|) but 
produce different data for the long-range spatial boson 
correlations. Here we limited ourselves to considering 
the most simple variant of generalizing Eq. (12), which 
makes it possible to investigate only the thermodynamics 
of a strongly interacting Bose gas. The interesting and 
important problem of the spectrum of the elementary ex- 
citations is thus beyond the scope of this paper and will 
be the subject of the future investigations. 



V. CONCLUSION 

Concluding let us take notice of the important points of 
this paper once more. It was demonstrated that thermo- 
dynamically consistent calculations based on Eqs. ( [ll| ) 
and (|l2|) conclusively result in Eqs. (|o|)-(||). There- 
fore, using the Bogoliubov model with the "dressed" in- 
teraction does not provide the satisfactory solution of 



the problem of the strong-coupling Bose gas. As it was 
shown, when investigating this subject, one should go be- 
yond the Bogoliubov scheme. To do this, we developed 
the approach reduced to the system of Eqs. (^3|), (|3S|), 
( ^ ) and (^) . These equations leading to the in- medium 
Lippmann-Schwinger equation (^), reproduce the famil- 
iar results (|4^), (|o|) and ( |5^ ) for the condensate deple- 
tion, chemical potential and mean energy but yield com- 
pletely different picture of the spatial boson correlations. 
This difference should manifest itself in the next orders 
of the density series expansions for the thermodynamic 
quantities and in the excitation spectrum as well. 
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